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In this article, we perform a careful analysis of the renormalization pro- 
cedure used in existing calculations to derive critical exponents for the KPZ- 
equation at 2-loop order. This analysis explains the discrepancies between 
the results of the different groups. The correct critical exponents in d = 2 + e 
dimensions at the crossover between weak- and strong-coupling regime are 
X = C(£ 3 ) and z = 2 + 0(e 3 ). No strong-coupling fixed point exists at 2-loop 
order. 
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I. INTRODUCTION 



The Kardar-Parisi-Zhang equation |T 
dh(x, t) 



A 



dt 



vV 2 h(x, t) + - (Vh(x, t)y + ri(x, t) 



T](x, t)r)(x', f) = 2D5 d (x - x')5{t - f) 



(1) 
(2) 



plays a central role as simplest field-theoretic model for non-linear growth and has extensively 
been studied during the last years. Whereas the situation is clear for space-dimension d = 1, 
the (physical more interesting) case of d > 2 can only be attacked by approximative methods 
or field-theoretic perturbative expansions. Using the latter, the fixed point structure of 
the renormalization group flow for d = 2 + e has been obtained. Two domains can be 
distinguished: For small effective coupling 



AD 



(3) 
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the renormalization group flow tends to in the long- wavelength limit, for large coupling to 
a strong coupling fixed point g = g sc . The crossover takes place at g = g co , which turns out 
to be of order e in an e-expansion and can therefore be studied perturbatively. 

Three such calculations have been performed, one by Sun and Plischke 0], another by 
Teodorovich |J, the third by Frey and Tauber [|J. The results are contradictory. We were 
able to trace back the difference to a problem in the application of the renormalization group 
procedure, which lead and |3| to incorrect results, see the discussion in sections fTT], f[V] 
and [V[ We will first derive some elementary and hopefully well-known properties of the 
renormalization group, which will help us to make the point clear. 



II. SOME ELEMENTARY PROPERTIES ABOUT RENORMALIZATION AND 
THE FREEDOM TO CHOOSE A RENORMALIZATION SCHEME 

In the following, we want to exploit the freedom one has to choose renormalization 
factors in the dimensional regularization scheme. As an example, take any massless theory 
and suppose that the derivative with respect to k 2 of the bare and renormalized 2-point 
functions (we denote this by a dot) at a given renormalization scale \x are related by 

f g> = Z(g R ) if . (4) 

The factor Z(g R ), which is a function of the renormalized coupling g R and s only, is intro- 
duced to define a finite renormalized two-point function T R \ Then, it is standard to obtain 
the anomalous dimension of the field, given by the renormalization group ^-function, see e.g. 
H, as variation of ln(Z) with respect to the renormalization scale /i, keeping bare quantities 
fixed 

1 d 

Cm) = - ^-qZ foZ(g R ) 

^ B 

= -l(3(g R )^-lnZ(g R ) . (5) 

In the second line, we have introduced the renormalization group /^-function, which is ob- 
tained from the variation of the renormalized coupling g R with respect to the renormalization 
scale fi: 



g R (6) 



B 



and where renormalized and bare coupling constants g R and g R are related by 

9b = li e Z~ x (g R )g R . (7) 
The critical point g* is the zero of the /3-function 

W) = • (8) 

The critical exponent (* is 
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c* = con . (9) 

All these definitions are standard || and are used (in slightly different notations) in @, || 
and 0.(0 uses a cut-off regularization, but there is no principal difference.) We now ask 
the question: Is the renormalization factor Z(g R ), which was introduced in equation (f|), 
unique? The answer is no. Z(g R ) can be multiplied by any finite factor Zf(g R ) and T R 
rests finite. Instead of Z(g R ) one could therefore use Z'(g R ), 

Z\g R ) = Z(g R ) x Z f {g R ) . (10) 

Let us verify that this is consistent and that the critical exponent (* is unchanged. The new 
C'-function is 

£{g R ) = -\^^z>{g R ) 

= ag R )-~(3(g R )-^lnZ f (g R ) . (11) 

As Zf(g R ) is finite for any g R , (3{g R )-^ \nZf(g R ) tends to for g R — > g* and we obtain the 
result 

cv) = con ■ (12) 

One can also prove that the freedom to choose g R is a simple finite reparametrization 

g' R = g R x Z f (jg R ) . (13) 

One has however to be careful to make this change of variables in every renormalization 
factor. 

(2) 

Let us now discuss the minimal subtraction scheme. Suppose that T B at the renormal- 
ization point /i satisfies up to 2-loop order, i.e. up to order g\ 

fg } = (l + + a g R + a l9R e + ... + 6_ 2 |f + + b g 2 R + hg 2 R e + (14) 

where and 6j are constants and where we have suppressed higher order terms in e. We 
now want to introduce renormalization factors. Of course, setting in (U) 

Z(g R ) = 1 + a^ 9 — + a g R + a x g R e + ... + + + b g 2 R + b ig 2 R e + ... (15) 

e e z e 

would render finite. (The dots stand for the same higher order terms as in (|I~4"D.) But 
this is not the simplest possible choice. The idea of the minimal subtraction scheme is now, 
to choose a Z-factor, which contains only the pole-terms in e. One is therefore tempted to 
replace Z by 

Z PM (9R) = 1 + a-i— + b-2?i + 6_i^ . (16) 

e e 2 e 
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The reader is invited to stop for a second to think about this choice before going on to read. 

Let us see whether Zpu correctly renormalizes the theory. Then, the ratio of Z and Zpu 
has to be finite. Expanding up to order g 2 R and dropping all terms of order e° and higher, 
we obtain: 

z (sr) _, aoa _ & + .... (17) 



This ratio is not finite and Zpm therefore not a correct renormalization factor! We will 
henceforth call this wrong prescription to find the Z-factor in the minimal subtraction scheme 
"pseudo minimal-subtraction scheme" in order to make clear that this very special mistake 
was done. 

The correct factor to be used in a minimal subtraction scheme would be 

Z M s(9r) = l + a^ + bj\ + (&_! - a a_!) ^ (18) 
e e z e 

as the reader may easily convince himself by the construction given above. 

Apparently, the term cannot be neglected at 2-loop order. Let us determine more 
generally, which terms may intervene in a n-loop calculation. First of all, the perturbative 
fixed point in an e-expansion is of order e 

g*~ £ . (19) 

(There are situations where g* is of order e 2 or higher. For those cases the following argument 
has to be modified.) Let us now determine in which order of e a term of order 

9 -f (20) 

in Z(g R ) will contribute to the (^-function at g R = g*. Taking care of the factor e, which 
comes from the variation of lnZ(g^) in the definition of ({go), equation fl5|), the answer is 
with the help of (|TJ) 

e (£!! = e i+*-J . (21) 

To calculate the critical exponent (* up to order e n , one therefore has to determine all 
contributions to Z which are proportional to g\ up to order e n ~ k ~ 1 . 

This means that to calculate up to order e 2 , one can drop terms from Z and retain only 

Z d (g R ) = 1 + + a g R + bj\ + bj^ . (22) 

e e z e 

Note again that a cannot be dropped. 

In the next sections we will apply these considerations to the 2-loop calculations of the 
KPZ-equation in the literature. 
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III. ANALYSIS OF THE WORK BY SUN AND PLISCHKE 



In this section we are going to analyse the 2-loop renormalization group calculations by 
Sun and Plischke ||, see also the comment in |J. From their Appendixes A and B it is clear 
that the finite part of the 1-loop contribution is dropped. They write: "In this Appendix 
we list all expressions of Feynman graphs . . . and their final results in 2 — e dimensions to 
order 0(e -1 )." As there is a finite contribution at 1-loop order, the discussion given above 
shows that their results are incorrect. We would have liked to incorporate the missing term 
in their calculations but as not all Z-factors are given explicitly we have to leave it to the 
authors. 



IV. ANALYSIS OF THE WORK BY TEODOROVICH 



Let us also look at the article by Teodorovich ||. Teodorovich is conscious that the finite 
part of the 1-loop diagrams enters as insertion of these diagrams into 2-loop diagrams and 
he explicitly calculates the finite contributions to the 1-loop diagrams in his equations (5.1) 
and (5.2). But then he also drops the finite contribution to which is still present in 
equation (5.1), in equation (5.3). We are therefore led to the same conclusion as in the case 
of Sun and Plischke: The results are incorrect. Let us however mention that Teodorovich 
does not agree with our criticism. He states to have used the BPHZ procedure and inserted 
1-loop renormalized terms in order to calculate the renormalization factors |]]. It is however 
impossible from his article to get the precise procedure enabling one to verify his claim. 



V. ANALYSIS OF THE WORK BY FREY AND TAUBER 



In this section we are going to analyse the 2-loop renormalization group calculations by 
Frey and Tauber j|, see also their comment in [0. 

They introduce the following two renormalization constants for the renormalization of 
D and z/, see j| equation (3.43) and (3.44) 



Z D 

z v 



e ' e 



1 + 



d-2 



9b 



+ {d 



d 2e 
g% d-2 g% 
] 2e d 2e 



(d-1) 



9% 
2e 2 



d 2e 



(23) 
(24) 



The notation is 



9b = if 9 b (25) 
e = d-2. (26) 

(Note that our cjb = go i n the notation of [|J. To simplify the considerations, we also have 
replaced the function F u (d) in (3.44) of |§ by its value at d = 2, see (A38) of gj. The 
difference is a term of order e° and can be neglected in 2-loop order.) As we will check 
later these Z-factors subtract all divergences correctly. Let us first explain the appearance 
of factors d and e. Frey and Tauber Q use a dimensional regularization prescription, where 
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divergences appear as poles in 1/e. This is the origin of the factors 1/e and 1/e 2 in (^) and 
(|||). On the other hand, there are integrals of the type 

J (2(kp) 2 - k 2 p 2 ) f(k 2 ) (27) 



which are identical to 

p2 ^r// /(fc2) - (28) 

Frey and Tauber call the factors (d — 2) "geometrical" because they result from scalar prod- 
ucts in the integrand. The renormalization scheme in Ref. @] is such that those geometrical 
terms are not e expanded in order to meet the condition that a fluctuation-dissipation the- 
orem is valid in d — 1 dimensions. This procedure results in the appearance of the first 
non-trivial term in Z v (see equation 



d — 2q B , . 

— — — . 29 

d e 

From the reasoning in |4| one could gain the impression that such a term can be dropped, 
if the distinction between factors d — 2 and e is not made. This would be incorrect, since 
- as we have explained in detail in section |I| - this term has to be retained in order to 
renormalize the theory. We would write fl29D as 



\ + 0(e)) 9b (30) 

and using our findings at the end of section |TI] conclude that the term of order e° has to be 
taken into account. (To be precise: In section |I| we had argued on the level of the Z-factors 
as functions of the renormalized coupling g R . As 9b = 9r + 0(,9r) this does not affect our 
reasoning for the one- loop term in (|29|) and (|30D.) 

This analysis shows that the distinction between factors d — 2 and e is not necessary in 
order to renormalize the theory above d = 2 dimensions. Especially, one may not misinter- 
pret |3| in the way that one could choose d fixed in the so called "geometrical" factors and 
then perform the limit e — ► elsewhere. This would result in a divergent term of order 1/e 
in the renormalization group functions, violating the very principle of renormalizability. 

The correct procedure to follow is the following: Expand the Z-factors in e only. Retain 
terms of order g B , §b/£, 9b/ £ an d 9b/ £2 - (Of course, there would be no harm in keeping 
higher order terms in e, but this is unnecessary.) Then calculate the renormalization group 
functions. 

We refrain from discussing the validity of a perturbative renormalization group approach 
below the lower critical dimension d = 2. Such a scheme would have to meet certain criteria, 
such as the validity of the fluctuation-dissipation theorem in d — 1. But, it is not clear to 
us how this procedure should be defined. 

In order to demonstrate the correctness of our procedure and in order to show that 
factors d — 2 and e have not to be distinguished, we use the Z factors of [Q in our scheme 
to establish the critical exponents. In the view of our discussion in section we first check 
that the finite term in the factor Z^, is correctly taken into account. This is trivially true 
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as normalizations in |4]] are chosen such that the only integral which has to be calculated in 
1-loop order gives exactly 1/e and no finite contribution. 

Replacing d = 2 + e and retaining only terms of order g B , c/b/e, 9b/ £ an d 9b/ £2 > the 
Z-f actors are: 



(31) 
(32) 
(33) 



Z d (9b) 


= 1 - 


9_B + 

e 


9l 

e 2 




Zu(9b) 


= 1 + 


1 

- 2 9B~ 


1 

" 2 


9% 

e 


z g {g B ) 


= 1 - 


f 3 

U+ 


I) 


9b 



- + -*)9b- 



e e A 



The last factor, Z s , was calculated via the relation 



(34) 



as stated in ||. It is now necessary to transform to renormalized quantities. Replacing g B 
in (|3T|) , (|32|) and ([33|) by (note that this is sufficient in 2- loop order) 



1 3 

: 2 



ffs = <?« + ( 7 + £ ) 9% + 0{g%) 



we obtain 



1 



9r 3 5f R 2 



2 e 



Zd{9r) 
Zu{9r) 

z g (g R ) 

The renormalization-group functions as defined in M are 

d 



l + - 2 9n 

« 3 g R 

1 Or . 

T e 



(35) 

(36) 
(37) 
(38) 



Cd(j9i 



d/j, 



\nZ D (g R ) 



B 



d 



P(9R)-^—foZ D (g R ) 



dg 



R 



(39) 



d 



\nZ v (g 



R 



B 







/3(gR)^—lnZ v (g R ) 



dg 



(40) 



(3(9i 



d_ 



9r 

B 



£ 9r 



1 -g R g^\nZ g (g R ) 

9r ( £ + (d(9r) -Ku(9r)) 



(41) 
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Using our results from above, we obtain the following 2-loop renormalizat ion-group functions 

3 o 



Cd(9r) 
P(9r) 



-Qr - 2^R 



e 



1 



9r 



2 

£9r 



9l 



1 + ^ ) 9r 



(42) 
(43) 
(44) 



which are identical to the results obtained in equations (3.62)-(3.64) of ||] after a full e 
expansion. There is one non-trivial (g* ^ 0) zero of the /^-function. It is located at 
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(45) 



Note that no fixed point of order 1 ( "non-perturbative fixed point" ) exists in contrast to 
and ||. This yields the critical exponents 



Cd 

c 



-e + 0(s 3 ) 



(46) 
(47) 



The more standard quantities, the roughness-exponent x an d the dynamical exponent z are 



X 

z 



0(s 3 ) 
2 + 0(e 3 ) 



In addition, the correction to scaling exponent 

d 



to 



dg 



-(3(9 



R, 



R 



9R=a 



is easily read off from ( f44|) 



UJ 



e + 0(e 2 



(48) 
(49) 



(50) 



(51) 



We have given all calculations in this explicit form to convince the reader that the standard 
^-expansion is sufficient in order to obtain the critical exponents given in ((48|) and (|49| ). 

For all fans of the minimal subtraction scheme, we leave it as an exercise that via a 
variable-transformation and finite renormalizations one can transform to the following Z- 
factors in the minimal-subtraction scheme 



Z d {9r) 
Zu(9r) 



9r 



+ 0(9'r) 



0(g R ) 



(52) 



Z 9 (9R) = 1 - f + 0{gl) 



Using these Z-factors, one of course obtains the same critical exponents to order e 2 . 

We also found some minor misprints in [|]]. For the readers conveniance we give a list 
here [Hi. 
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VI. CONCLUSIONS 



In this article we have analyzed the 2-loop renormalization group calculations by Sun 
and Plischke 0, by Teodorovich and by Frey and Tauber J|. Sun and Plischke and 
Teodorovich || don't correctly take into account the finite part of the 1-loop diagrams. 
Therefore, their results are incorrect. 

The renormalization scheme by Frey and Tauber [Q leads to the correct result. From 
the reasoning given in their paper one could, however, gain the wrong impression that a 
distinction between socalled geometrical factors and \ je poles is necessary in order to get a 
valid renormalization of the KPZ equation. We have shown here, that such a distinction is in 
fact not necessary above the lower critical dimension d = 2. This conclusion is supplemented 
through work by Lassig, who showed [ID|, using the mapping to directed polymers, that the 
predicted critical exponents 

C = o 

z = 2 (53) 

are correct to all orders in e. This result had already been obtained in a different context 
in 
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On the technical level, we have shown that the dimensional regularization scheme works 
perfectly well and has not to be modified in order to calculate the critical exponents for the 
KPZ-equation to order e 2 above the lower critical dimension. The analysis of the critical 
behavior below d — 2 is more subtle and it is presently not clear whether there is a valid 
concept for the analysis of the strong coupling behavior. 
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